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A Study of the Afterglows of Four GRBs: Constraining the 
Explosion and Fireball Model 
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ABSTRACT 

We employ a fireball model of the gamma-ray burst explosion to constrain 
intrinsic and environmental parameters of four events with good broadband af- 
terglow data; GRB 970508, GRB 980329, GRB 980703, and GRB 000926. Using 
the standard assumptions of constant circumburst density and no evolution of 
the fraction of the explosion energy in the post-shock magnetic field, we investi- 
gate the uniformity of the derived explosion and shock physics parameters among 
these events. We find a variety of parameters: densities that range from those 
of the ISM to diffuse clouds, energies comparable to the total GRB 7-ray en- 
ergy, collimations from near-isotropy to 0.04 radians, substantial electron energy 
fractions of 10-30% with energy distribution indices of 2.1-2.9 and magnetic en- 
ergy fractions from 0.2-25%. We also investigate the level to which the data 
constrain the standard model assumptions, such as the magnetic field evolution, 
and the allowed density profiles of the medium. Fits generally improve slightly 
with an increasing magnetic energy fraction e#. Good fits can be produced with 
magnetic energy accumulating or decaying with the shock strength over the af- 
terglow as tB oc 7^; —2 < x < +1 . The data are not very sensitive to increasing 
density profiles, allowing good fits even with density oc r 10 . Some parameter 
values change by up to an order of magnitude under such altered assumptions; 
the parameters of even good fits cannot be taken at face value. The data are 
sensitive to decreasing densities; r~ 2 profiles may produce reasonable fits, but 
steeper profiles, even r~ 2 - 5 , will not fit the data. 
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1. Introduction 

Afterglow emission, spanning from radio to X-ray, and observable from minutes to weeks 
after the event, appears to be a ubiquitous feature of long-duration gamma-ray burst explo- 
sions. Afterglows result from radiation emitted by electrons accelerated in a relativistic shock 
produced by the explosion of a progenitor (e.g., Paczyhski & Rhoads 1993; Katz 1994). As 
the shock propagates outward into the surrounding medium, the resulting broadband syn- 
chrotron radiation evolves in a manner dependent on a number of fundamental characteristics 
of the explosion, as well as on the details of the shock evolution and the density profile of 
the medium it expands into (for a review of the theory, see Meszaros 2002). 

The physics of relativistic shocks, the mechanism of collisionless shock acceleration, and 
the means by which magnetic fields can be amplified to levels needed to produce the ob- 
served synchrotron emission are poorly understood. Basic understanding of these elements 
will eventually result from improved magnetohydrodynamic simulations, which may ulti- 
mately prescribe such factors as the appropriate electron energy distribution, what fraction 
of shock energy goes into magnetic field, and how this evolves with the shock expansion. 
Additional uncertainties arise from the fact that the outflow geometry, and the structure of 
the circumburst medium are unobservable, and could potentially be complex. 

In spite of the lack of detailed understanding, the basic features of the GRB afterglow 
can be described by relatively simple theoretical models for the outflow and shock (e.g. 
with relativistic flow as in Blandford & McKee 1976). Simple assumptions about the shock 
microphysics and geometry appear to fit the basic features of the afterglow in a number of 
cases (e.g., Wijers & Galama 1999; Chevalier & Li 2000; Harrison et at. 2001; Panaitescu 
& Kumar 2001b, a, 2002). In addition, many events appear consistent with the expansion of 
the shock into a medium with either a constant or a simple r~ 2 powerlaw in density. 

Ideally, high-quality broadband afterglow observations, interpreted in the context of 
basic theoretical models, could be used to constrain the explosion parameters, geometry, 
and the structure of the surrounding medium. In addition, if the data are of sufficiently 
high quality, it should be possible to test the validity of the basic model assumptions; or at 
least the range over which they may be varied and still describe the data. Given a relatively 
large number of model parameters, however, only a few data sets are of sufficient quality to 
provide interesting constraints. 

In this paper, we examine four high-quality broadband afterglow data sets, fitting them 
to a basic afterglow model. We investigate the similarity of the derived explosion, shock, 
and environmental parameters. In addition, we investigate selected model assumptions and 
the range over which they can be varied. In particular, we consider whether the fraction 
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of energy in the post-shock magnetic field can evolve as some power of the bulk Lorentz 
factor, and the range of possible matter density gradients. Finally, we consider which future 
observations provide the most promise for better constraining both the model and its physical 
parameters. 

2. The Afterglow Framework 

We adopt a standard fireball scenario for the GRB afterglow, where a relativistic shock 
with Lorentz factor 7 expands into the circumburst medium (CBM). Blandford & McKee 
(1976) describes the flow. The shock may be spherical or initially confined to a cone (i.e., 
be jet-like). The afterglow flux arises from the radiation (synchrotron and possibly also 
inverse Compton) emitted by relativistic electrons accelerated in the shock. To describe its 
evolution we model the shock dynamics, adopting assumptions for the density profile of the 
CBM, the shock geometry, and the shock microphysics including the distribution of electron 
energies, as described below. We also account for the effects of the medium through which 
the radiation passes en route to the observer. 

2.1. The Basic Afterglow Model 

For the dynamics, we make simple assumptions about the event's geometry and envi- 
ronment. We consider a constant circumburst matter density n, the "ISM-like" well 
as the possibility the CBM could be dominated by a wind outflow from the progenitor. As a 
constant mass-loss rate and wind speed gives an r~ 2 profile, we consider this the "Wind-like" 
case. Into either CBM form we allow for isotropy or a simple collimation of the ejecta, a 
top-hat distribution in solid angle with half-opening angle 0. The shock behaves as though 
isotropic until it slows down sufficiently to expand in its rest frame (Rhoads 1999; Sari, Piran 
& Halpern 1999). We calculate the time tj e t at which collimation is evident as 6 = I/7 (7 
the shock Lorentz factor; this assumes the observer is nearly along the line of sight). We 
calculate the expected radiative losses, which can modify the shock dynamics (Sari 1997), in 
the manner detailed further below. 

For the microphysics governing the emission it is standard to assume several things. 
First, that the shock imparts a constant fraction of its energy (e e ) to the swept-up electrons, 
and a constant fraction (e#) goes into amplifying magnetic fields. These fractions are capped 
at 100%, where all the shock's energy would be in one of these. A reasonable limit could 
have been the ejecta-e e -es equipartition value of 33%, but the model's calculations are not 
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perfectly known. The equations used have some uncertainties; we allow each of e e and 
€b to vary up to the limit of taking all the shock energy. It is also assumed that the 
electrons are accelerated into a simple powerlaw distribution of energies above a minimum 
value (P(7 e )oc 7~ p , 7 e > j m ), with a constant index p. 

These are our basic assumptions for the fireball model. They will be tested in further 
sections. 



2.2. The Emission calculation 

Given this electron spectrum, we get a broken powerlaw radiation spectrum, with three 
break frequencies that evolve in a manner that depends upon the dynamics (Sari, Piran & 
Narayan 1998): the injection break u rn for the minimal energy of the radiating electrons, 
the cooling break v c corresponding to energies at which radiative losses over the shock's 
lifetime are significant and the self-absorption break v a where the spectrum becomes optically 
thick at low frequencies. If the minimum energy electrons emit at a peak frequency above 
the self-absorbed regime, the spectrum below v a is oc u 2 , ressembling a blackbody with 
effective temperature corresponding to that minimum energy; if electrons are emitting within 
the optically thick regime, the effective temperature is a function of frequency and the 
spectrum is oc z/ 5 / 2 (Rybicki & Lightman 1979). We use a smooth shape with factors of (1 
+ (v / VbreakY 1 " 132 )" 1 (Pi, 2 are the indices before and after the break) for v m and u c , and for 
u a we use the physically motivated prescription of Granot, Piran & Sari (1999a). 

For the model with basic assumptions we employ the previously published fireball model 
calculations for the spectral breaks. The isotropic case is taken from Sari, Piran & Narayan 
(1998) 's relativistic near-adiabatic case, with normalizations from Granot & Sari (2002) 
(their canonical order, Table 2) for a Wind-like CBM, and from Sari, Piran & Narayan 
(1998); Granot, Piran & Sari (1999a,b) for an ISM-like CBM. These normalizations account 
for the post-shock spatial distribution of electrons and emission from equal arrival-time 
surfaces. The resulting equations are: 

F™ ax = 1.6(z + l)D 2 2 e° B ^_ 2 E 52 n - 5 mJy 

= 4.2xl0 8 (z + l)- 1 //(p)e e - 1 e^_ 2 ^ 5 °2 2 n a6 Hz 

v m 3.3xl0 14 (^ + l) a5 4 5 - 2 ^^5° 2 5 ^ L5 Hz [i) 

v c = 6.3xl0 15 {z + l)~ - 5 e B ^ 2 E^- b n- l t d °- Fj Hz 

For the ISM-like CBM, where f T (p) = ( (p + 2)(p- l)/(3p + 2) ) a6 and e e = e e (p- 2)/(p- 1). 
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For the Wind-like CBM: 

F™ ax = 7.7 (p + 0.12) (z + iy^D^e^E^A,^ - 5 mJy 

"« = 3.3 x 10 9 (z + I)" - 4 f w (p) 6- 1 e B % A\- 2 ff* Hz 

i/ m = 4.0 x 10 u (p-0.Q9)(z + l) - 5 e B 5 _ 2 e 2 e E° 52 5 tf- 5 Hz l J 

z/ c = 4.4 x 10 13 (3.45 - p)e 0A5p (z + l)" 1 - 5 'e B ^ 2 E° 52 5 A~ 2 t° d 5 Hz 

where fwip) — ( (p ~~ l)/(3p + 2) ) a6 . The units are: D 2 % = luminosity distance in 10 28 cm, 
€b-2 = 6b in %, n = density in cm -3 , A* = density scaling for the Wind-like profile so that 
p = 5 x 10 11 A*r~ 2 g cm -1 (r in cm; a standard reference for a mass loss of 10 -5 M Q yr _1 
at a wind speed of 1000 km s _1 ), E 52 = isotropic-equivalent energy in units of 10 52 ergs, td 
= observed time post-burst in days. The electron energy partition e e is given as a fraction. 

The post-tj e t evolution with laterally expanding ejecta uses the time-dependences of 
Sari, Piran & Halpern (1999) (eqs. 2-5), and is connected to the pre-jet behaviour without 
smoothing (i.e., a sharp jet break). 

Using the relativistic equations for the shock energy, E = M7 2 , we calculate the time 
at which 7 = 1 as the nonrelativistic transition, t nr . This would be equivalent to using the 
Blandford & McKee (1976) approximation for the energy, E = M^y 2 (3 2 , and defining the 
nonrelativistic transition condition as 7/9 = 1. We again employ a sharp transition to the 
post-t nr behaviour. 

We adjust this evolving synchrotron spectrum with self-consistent corrections to the 
cooling rate, based upon the parameters already enumerated. We calculate the effects of 
synchrotron photon upscatters (inverse Compton scatters, IC) off the shocked electrons 
as in Sari & Esin (2001). When IC dominates the cooling we adjust the cooling break 
v c (by (1 + Y)~ 2 ~ Y~ 2 when Y > 1; Y the IC to synchrotron luminosity ratio R3 
(max[(7 c /7 m ) 2_p ,l]e e /e5) 1 / 2 ). The upscattered photons can produce a high-frequency sec- 
ondary source of flux; for it we adopt the spectral breaks as in Sari & Esin (2001) and use 
our synchrotron spectral shape (a good approximation except for an ignored logarithmic 
correction between v 1 ^ and v 1 ^ and the slope for v < z/^ c , where IC never dominates the 
data). 

We also treat radiative corrections self-consistently, from the model parameters. In- 
stantaneously, we treat the shock as adiabatic. The energy is calculated for each time from 
the solution of dE/E (Cohen, Piran & Sari 1998). It depends upon the ratio u c /u m in slow 
cooling [y c > u m ) as the losses quench out and we use the synchrotron-only rate of change in 
z/ c , allowing a simple analytic solution for E{t). As this is when E(t) changes are becoming 
unimportant, the approximation has little effect. We scale the shock's energy to the value 
at the change from fast- to slow-cooling regimes [y c = u m , fairly early). 
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This spectrum is modified on the way to the observer, and we must account for extinc- 
tion, interstellar scintillation (ISS) and the host's flux. ISS in the Galaxy distorts the flux 
at low (radio) frequencies; we estimate its fractional flux variations as outlined in Walker 
(1998) for point sources, using the map of scattering strengths by Taylor & Cordes (1993), 
and scalings for extended sources as explained in Narayan (1992). The growing angular size 
quenches the scintillations; we assume initial parameters in the model to get angular size as 
a function of time for the ISS model. These are not iterated, but used as an additional un- 
certainty in the model flux, added in quadrature to the data's uncertainties when estimating 

x 2 . 

For extinction, we first deredden for Galactic effects using Schlegel, Finkbeiner & Davis 
(1998) 's method for estimating A A in the optical. As most GRBs are at high Galactic lati- 
tudes, this term is small. Then, as the extinction curves in high-z galaxies are unknown, we 
use known local extinction curves redshifted to the host restframe to fit for host contribu- 
tions. We use Weingartner & Draine (2001)'s Large Magellanic Cloud (LMC) in general, but 
for the GRB 980329 event we used the Small Magellanic Cloud Bar's law due to evidence 
that the host's extinction law is steep. We do not subtract host fluxes as the decay slope is 
sensitive to small differences in the host value; host fluxes are added to the model in the fit. 
In the optical and near-IR, when there is evidence for a host component, we fit a value for 
the particular band(s) involved. As submillimeter hosts have been detected in other bursts 
(Hanlon et al. 2000; Berger, Kulkarni & Frail 2001; Berger et al. 2001; Frail et al. 2002; 
Berger et al. 2003), we allowed for this possibility (scaled to 350 GHz as z/ 3 ), but none was 
required for the four events under consideration and this component was not included in the 
best fits. Finally, we see in 2 datasets evidence for underlying flux in the radio. 980703 has 
sufficient data that a spectral index (-0.32, Berger, Kulkarni & Frail 2001) was determined 
and we fit the radio host flux scaled to 1.43 GHz (where it is brightest). For 980329 we use 
a canonical spectral index of z/~°- 8 , again scaled to 1.43 GHz. 

Two other groups have performed full analyses similar to ours, as described in Chevalier 
& Li (2000); Li & Chevalier (2001) and Panaitescu & Kumar (2001b,a, 2002). There are 
some differences in each group's approach to the fireball model. Chevalier & Li develop an 
r~ 2 CBM model, including analytic solutions to the hydrodynamics, but they do a numerical 
solution to get the full smooth spectral shape. They assume an extra break in the electron 
energy distribution but do not, however, account for energy losses, IC cooling or ISS effects. 
Panaitescu & Kumar numerically solve the dynamics (whereas we use analytic asymptotic 
forms), numerically calculate the spectrum from equal arrival-time surfaces (whereas we use 
equations adjusted to account for equal arrival-time surfaces and the electron distribution 
behind the shock, and a simple smoothing) and calculate IC scatters and energy losses 
directly from the model's radiated spectrum (whereas we use the theoretically expected levels 
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from the parameters). They use the simplest electron energy distribution when it provides 
a good solution but also allow for two (rather than one) indices in the injected spectrum if 
needed. They pre-account for host effects by subtracting fluxes and fixing extinction levels, 
instead of including these in the fits. 

3. The Sample of Four GRBs 

We selected four events with rich datasets (radio, optical-NIR and X-ray at a minimum) 
for which the fireball model with simple assumptions provides a good description of the 
lightcurves. These are GRB 970508, detected by the BeppoSAX GRBM instrument, with its 
X-ray afterglow found by the BeppoSAX WFC (Costa et al. 1997), GRB 980329, its afterglow 
detected in the X-rays by In 't Zand et al. (1998), GRB 980703, detected by BATSE (Kippen 
et al. 1998) with its afterglow identified by the Rossi X-Ray Timing Explorer (Levine, 
Morgan & Muno 1998) and GRB 000926, located by the IPN (Hurley et al. 2000; Hurley 
2000), with its afterglow first identified in the optical (Gorosabel et al. 2000; Dall et al. 
2000). 

We compiled and described the afterglow data for the 980329, 980703 and 000926 events 
in previous papers. The summaries presented here chiefly point the reader to the relevant 
papers. As we have not compiled the 970508 data elsewhere, we provide full details here. 

GRB 970508 was the second burst with a detected afterglow, and it has a rich broadband 
dataset. We use the X-ray fluxes of Piro et al. (1998), converting to flux densities and 
frequency assuming v~ 1A over their quoted band. For the optical/NIR, we use BVI data from 
Sokolov et al. (1998) and Zharikov, Sokolov & Baryshev (1998), R from the table compiled 
in Garcia et al. (1998), and K s from Chary et al. (1998). We converted magnitudes to flux 
densities with the information in Tokunaga (2000) for K s and Bessell (1979) for BVRI. We 
included the 12- jum mid-IR measurement by Hanlon et al. (1999). We used Bremer et al. 
(1998) in the submillimeter (86 & 232 GHz), and the 15 GHz data of Taylor et al. (1997) 
and the 1.43, 4.86 and 8.46 GHz dataset of Frail, Waxman & Kulkarni (2000) in the radio. 

GRB 970508's afterglow was unusual - a sudden brightening in (at least) the optical 
occurred at ~ 1 day, followed by the usual powerlaw decays at t > 2 days, and the dataset 
shows an intrinsic "scatter" even at later times in at least the optical and radio. No simple 
model can account for the rise. Various possible explanations for a sudden rise in afterglow 
flux include a refreshed shock from slower material in the relativistic flow catching up to the 
slowing forward shock (Panaitescu, Meszaros & Rees 1998; Sari & Meszaros 2000; Kumar 
& Piran 2000), a sudden jump in the encountered density (Wang & Loeb 2000; Dai & Lu 
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2002; Nakar, Piran & Granot 2003), or even the "patchy shell" model where there are hot 
and cold spots in the shock and as the shell slows a hot spot may come into the observable 
area (Nakar, Piran & Granot 2003). These models can provide a single rise that will not 
affect the emission substantially later. To avoid such complications we used only the data 
at t > 2 days, after the rise, for the modelling dataset. 

The 980329 event has data in the X-ray, optical, submillimeter and radio regimes. Initial 
searches in the optical were unsuccessful, and Taylor et al. (1998) first identified the afterglow 
in the radio. Early epoch optical data where the afterglow dominates over host emission is 
therefore somewhat sparse. The excess t < 3 days emission at 8.46 GHz is thought to be due 
to a reverse shock (Sari & Piran 1999b) and therefore we excluded it from our analysis. Yost 
et al. (2002) gives full details concerning the dataset. No redshift has been determinable 
for this event (the only such case in our sample) due to faint host emission and a lack of 
prominent emission lines; we adopt z=2 here, roughly the middle of the potential redshift 
range. 

The 980703 afterglow data includes radio, optical/NIR (with good frequency coverage: 
BVRIJHK) and X-ray. The utility of the opt-NIR data is limited to a few days postburst 
by a bright host (which is also visible in the late radio data). Frail et al. (2003) give the 
details. 

The 000926 afterglow data includes radio data from 1.43 to 98 GHz, ground-based 
and HST optical and NIR observations and ChandraX-r&y observations. The first 8.46 GHz 
point (~ 1.2 days postburst) was excluded from the analysis as abnormally high and possibly 
associated with reverse shock emission. We use all available BVRHK data and divide the 
X-ray into a soft and a hard band. This is further described in Harrison et al. (2001). 

4. Fits to the Basic Model 

970508 shows a "scatter" in its lightcurves in all frequency regimes and timescales. It 
is not clear how much of the scatter may be due to real physical effects (such as dumpi- 
ness in the circumburst medium), and how much could be attributed to cross-calibration 
uncertainties. While we performed our fits as in §2, the best fit, passing through the data, 
has x 2 /DOF w 2.5. (We do not attempt an analysis with increased error bars to reduce x 2 
to the DOF, as this would require much larger uncertainties.) We have adopted this fit, a 
nearly-isotropic solution (tj et = 183 days) with an ISM-like CBM of somewhat low density 
(0.2 cm -3 ) and e e and e# near equipartition, as our best solution. 
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Yost et al. (2002) fit 980329 as detailed in §2. Our best model is collimated, with a 
significant density (n ~ 20 cm" 3 ) and host extinction. We found no need for a sub millimeter 
host component, and that a radio host improves the fit with marginal significance. The best 
fit was not a unique solution; Yost et al. (2002) noted an isotropic, extremely radiative {all 
the shock's energy in radiating electrons) solution fit as well. We reject that solution on the 
basis of its unphysical parameters. 

Frail et al. (2003) gives results of 980703's fits as in §2. The best fit is a collimated jet 
into a constant density CBM, its jet break hidden in the optical by the host's dominance. It 
was not a unique solution; a collimated jet in an r~ 2 CBM fit the data equally well. That 
fit required extreme parameters: 70% of the shock's energy is imparted to the electrons and 
the energy during the afterglow is ^ 1/10 that of the prompt 7-ray emission, requiring an 
extremely high 7-ray conversion efficiency We reject the r~ 2 fit. 

The model we employ (§2) for the 000926 fits has changed since the work done in 
Harrison et al. (2001). Before, estimated host fluxes were subtracted instead of fitted for. As 
there was a galaxy "arc" nearby contaminating the ground-based data, the host components 
fitted for are merely the persistent underlying flux common to both ground- and space-based 
observations; the "arc" flux removal remains as before. Moreover, the model now employs 
radiative corrections to the energy, which were not yet incorporated in the previous analysis. 
We refit the dataset; the best fit from the improved model is in Table 1. It is still a high-IC 
ISM model with a substantial density, but the x 2 has increased. The previous fit had a 
nonnegligeable e e and no radiative corrections; with radiative corrections correlated with e e , 
these are incompatible, e e drops and the fit worsens. In particular, while the X-ray still 
requires flux above the synchrotron's level, the parameter changes cause the IC flux estimate 
to no longer match the data as well. 



4.1. Comparisons of the Four Bursts 

Table 1 presents the best fits using the basic assumptions described in §2.1 (also see 
Figure 4). The table includes statistical 68.3% confidence intervals for the parameters, 
calculated from (non-Gaussian) distribution histograms generated by over 1000 Monte Carlo 
bootstraps. The results show a great deal of diversity in the values, similar to other efforts 
(e.g. Panaitescu & Kumar 2001b). The fits can be seen in Figures 1, 2, 3. 

Neither the variation in nor the values of the energy/geometry and environmental pa- 
rameters is surprising. Variation is seen in the energies of the prompt GRB emission (Frail 
et al. 2001). Moreover, if GRBs are related to the deaths of massive stars (see the review 
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by Meszaros 2002, and references therein) then we would not expect such parameters to be 
identical due to variations in progenitor mass, angular momentum and environment. 

The densities are comparable to the Milky Way's ISM density at the low end, and to the 
density of diffuse clouds for the three ~ 20 cm~ 3 , typical of other efforts with significant radio 
data (which constrains the synchrotron self-absorption) (e.g., Panaitescu & Kumar 2001a). 
The values are not inconsistent with the core-collapse GRB progenitor hypothesis, despite the 
expectation of a stellar birth environment, and molecular cloud cores having n 10 4 cm~ 3 . 
Diffuse clouds are also associated with star-forming regions, and such densities are found in 
the interclump medium that dominates molecular cloud volumes in the Galaxy (Chevalier 
1999, and its references). Moreover, massive stars modify their environment making direct 
interactions with the dense cloud cores unlikely (Chevalier 1999). 

The total kinetic energies inferred in the fits (10 51 -10 52 ergs) are roughly comparable, 
though at the high end, with the range of total 7-ray energies expected from the studies by 
(Frail et al. 2001). On a case-by-case basis, we check if they are reasonable by comparison 
with the 7-ray energies, employing the isotropic-equivalent k-corrected 20-2000 keV fluences 
of Bloom, Frail & Sari (2001). The ratio of the isotropic-equivalent kinetic (at 1 day post- 
burst) and 7 (over the GRB) energies varies from 0.5 to 5. However, the modelled energy 
corrections would give fireball energies 3 to 10 times higher at the end of the burst; the 
implied 7-ray efficiency would then be 3-10% in three cases, and 40% for GRB 000926. The 
lower efficiency values are quite compatible with theoretical expectations; 000926's is near 
the high end of the predictions (see Beloborodov 2000; Spada, Panaitescu & Meszaros 2000; 
Kobayashi & Sari 2001). 

Thus the environmental and energy parameters are quite reasonable in the present 
framework. What is somewhat unexpected is the lack of universality in the microphysical 
parameters. Values of e e are fairly uniform, only varying by a factor of 2, but the values 
of €3 vary by a factor of 100. The index p is expected to be in the range of 2.2-2.3 (see, 
e.g., Kirk et al. 2000; Achterberg et al. 2001) but its values span from 2.1 to 2.9. While 
relativistic shocks are not well-understood (neither fully modelled from first principles, nor 
measured in the lab), the physics occurring at a shock boundary should only be a function of 
the shock strength, or equivalently here for a relativistic shock, its Lorentz factor. A spread 
by two decades in e# would therefore be unexpected if everything is properly accounted for 
in the model. We would expect the microphysical parameters to be near-universal. 

If the spread in the microphysical parameters does indeed exist then highly relativistic 
shocks behave nonintuitively. (However, the apparent spread could be due to the overall 
uncertainty in the model, a parameter such as the viewing angle that is not fully accounted 
for or nonuniqueness of the model fits.) We have already noted that there are at least 
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two nonunique fits: GRB 980329 with a highly radiative fit, and GRB 980703 with an r~ 2 
density profile fit disfavoured due to extreme parameters. This leaves open the question 
as to whether there may be equally good fits with reasonable parameters (perhaps under 
other model assumptions). It is therefore important to check if the model, with the very 
simple assumptions we have employed, is constrained only to these fits by the data. The 
assumptions could be too simple, and the following section explores fits with other model 
assumptions in order to constrain some interesting aspects of the model uncertainty. 



5. Constraints on Deviations From the Basic Model 

Above, we found good fits even with simple assumptions. But without quantifying the 
uncertainty in model assumptions, it is not clear if this says anything about the simplicity 
of the underlying physical processes or the environment. We do not know if fits are possible 
with a variety of underlying assumptions, and whether that would substantially change the 
parameters inferred. We now investigate the range of allowed modifications on a limited set 
of fireball model assumptions, and the effect of modified assumptions on the parameters in 
good fits. 

First, the magnetic field amplification (assumed to be from instabilities) is not well 
understood - we do not know from first principles whether the fields pile up, decay away or 
reach a steady state. There is also a spread in the magnetic energy fraction e# (by a factor 
of more than 100). So the behaviour of the magnetic fields is clearly one of the important 
questions surrounding the microphysics of relativistic shocks. We want to learn whether it is 
possible for the field energy to evolve with the strength of the shock in the afterglow phase, 
and if the data truly constrains the interesting diversity seen in the best fit values. To get 
some constraint on this §5.1 explores whether an evolving €b could fit the data. 

Second, we know there is at least one case where more than one density profile can 
produce a good fit (§4). There is also a growing body of evidence that GRBs are associated 
with star-forming regions and may be the result of the collapse of massive stars (see Meszaros 
2002, and references therein). If correct, the shock should be propagating into a medium 
enriched by the stellar wind of the progenitor, and the density profile should not be constant. 
However, the best fits are for a constant-density medium. We want a constraint on how 
sensitive the fits are to density profile, and so we investigate a wider variety of density 
profiles (p oc r s ) than the two previously considered, in §5.2. 
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5.1. Magnetic Energy Fraction as Function of ^shock 

The strength of the magnetic fields implied in the fireball model are far stronger than the 
levels expected from the strength of the Galactic fields. The relativistic shock is expected 
to amplify the nearby fields (e.g. Medvedev & Loeb 1999). This would be a nonlinear 
process acting on small instabilities; the resulting fields would depend on the amplification 
mechanism and not on the initial values. This aspect of relativistic shock physics is not fully 
understood (see the review by Kirk & Duffy 1999, and references therein). GRB afterglow 
emission can be used to investigate this process; Rossi & Rees (2003) show that spectra will 
differ if the magnetic field is confined to a small length behind the shock. The importance 
of magnetic fields in GRBs is further demonstrated by Coburn & Boggs (2003) 's recent 
detection of linear polarization at the theoretical maximum in the prompt emission of GRB 
021206. This would require a uniform magnetic field across the 7-ray emission region beamed 
to the observer, potentially implying that magnetic fields dominate the dynamics (Lyutikov, 
Pariev & Blandford 2003). 

The assumption that the magnetic energy fraction e# imparted by the shock amplifica- 
tion is constant is the simplest assumption; shock physics could depend on the Mach number 
(equivalently the Lorentz factor if relativistic). Therefore it is possible that the shock would 
have different efficiencies with time, and the fitted number is simply some averaging over 
where its modelled effect on the data is strongest. It would also be possible to get effectively 
an increase or decrease in e# if the pileup of magnetic fields behind the shock in the place 
where the electrons are accelerated and radiate does not reach a steady state - if the fields 
remain and grow it would increase and if they diffuse away faster than they are replenished 
it would decrease. (It is also therefore possible that the pileup could be position-dependent, 
even if there is a universal behaviour for a collimated flow, and some disparities could be 
due to the unaccounted-for differences in viewing angle.) 

It is important to constrain whether the data allows a variable ee to see if such expla- 
nations are viable, and its spread could be due to model uncertainty. If the e B disparities 
could be resolved with a universal nonconstant e B behaviour that would be of significant in- 
terest. Moreover, if e^'s behaviour is nonconstant, the values fitted with the simple constant 
assumption could be different from the real magnetic fraction; the level of this effect should 
be checked. And if only a growing e B (pileup of field energy as the shock strength drops) or 
a falling e B (drop in the field energy as the shock strength drops) are allowed by the data, 
this would be a clue to the unknown details of relativistic shock's field amplification. To see 
what the data permit, we studied a simple parametrization oi e B with the shock's 7. 

We took the equations for the model as presented in §2 and allowed the fixed magnetic 
energy fraction to vary smoothly. A value for the parameter is calculated at each time and 
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inserted into the equations governing the spectrum at that time. Under assumptions where 
€b grows we do not permit it to get infinitely large, capping it at 100%. Any change in 
the magnetic energy fraction is expected to come from the evolution of the shock, and thus 
the simplest physically motivated form is to tie it to the shock strength, as expressed by its 
Lorentz factor 7. For simplicity, we considered e# oc r y x , where the basic model (§2.1) has 
x = 0. 

This affects the results in two ways, directly and indirectly. There is an explicit depen- 
dence upon €b in the spectral breaks (Sari, Piran & Narayan 1998) whose evolution changes. 
Indirectly, as the cooling frequency affects the energy losses as described in §2.2, we also 
recalculate E(t), which affects the dynamics. 4 

The resulting equations follow [y a is for the order v a < u m < u c . When u c < u m , 
multiply by (u m /u c ) 1/2 ). 









t < tj e t 


tjet *^ t <C t nr 




Fmax OC 


£ n 0.5 e e 0.5 




^O-O.l&r 


£-l-0.25x 


t 0.6 


v a OC 


£0.2 n 0.6 e -l e 0.2 




^0-0.082; 


^-0.2-0.1x 


t 1.2 


v m oc 


E°- 5 n°e 2 e e° B 5 


t 


-l.5-0.19x 


£-2-0.25:r 


r 3 


v c oc 


E-o-Sn- 1 ^ 1 - 5 


t 


-0.5+0. 56x 


^0+0. 75x 


r 0.2 



We examined model assumptions of several indices x, x — 2, 1, 0, -1, -2, -3. We did not 
free-fit the index, considering instead integer values. Changes are not highly sensitive to the 
index; for a typical time span of a factor of 100, 7 changes by ^1/6-1/10, and the change 
in e B for x — 1 alters the magnetic field by only a factor «2.4. Table 2 gives the expected 
behaviour of the synchrotron flux for x ^ 0; these changes are not large for 5x <C 1. 

We thoroughly searched for the best fit models at the values of x considered. Our 
methods included stepping x in very small steps (between 0.05 and 0.2, depending upon the 
ease with which the fitting would adjust for a dataset), allowing small, smooth changes. We 



4 Note that for x < —4/3 synchrotron-only theory predicts a decreasing v c jv m requiring us to change 
the function E{t). The transition to fast cooling, where all injected electrons can radiate a significant part 
of their energy, would only permit energy losses to quench out once es stops changing - as late as the 
nonrelativistic transition. In practice this is not the case in models of interest with es rising over the data 
range; e e is large enough that e e > e_e- IC cooling dominates so v c w v p) nch x es/e e (Sari & Esin 2001). 
With moderate energy losses at even x = — 3 the ratio rises at a low rate. The radiative corrections fall 
between no quenching of the energy losses and the quenching rate of x — 0. These arc not very different; 
typically the energy difference between these out to ~ 200 days is ~ 2 x , a small effect on the final data with 
low S/N. We use the synchrotron v c jv m ratio to calculate quenching above x = —4/3 and check both no 
quenching of E{t) and full x = quenching for x < —4/3. These give substantively the same results. 
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tried larger steps 5x (~ 0.1 where 0.05 would work smoothly, ~ 0.25 or 1 where 0.2 would 
work for small steps), forcing the gradient search algorithm to look farther afield for a good 
fit. We also fit from grids of selected parameter starting points at a particular (integer) x. 
These grids included values comparable to those at the best fit for the x next-nearest to 0, 
the basic model's values and typically went up & down by a factor ~ 10. 

Best fits for various x are in Table 3. We find that x < fits as well or better than the 
basic model out to at least x — — 1; the best oc 7 _1 fits including uncertainties are shown 
in Table 4 and Figure 6. The improvement is seen in 3 of the 4 cases, while the fit with the 
x = — 1 model assumption to 970508 is only 1% worse in the total broadband x 2 than the 
basic, constant-es fit. The effect of the x — — 1 assumption on the fit to optical and X-ray 
data is minor. It takes a long time baseline, as in the radio, to see the spectral evolution 
differences; moreover the optical and X-ray are above the cooling frequency and thus have 
a lesser flux dependence upon e# (Table 2). Assuming es oc 7 _1 improves the radio fits - 
the magnetic energy grows as the shock slows, which allows the peak to rise, flattening the 
late time decay, as shown in Figure 5. We note there appears to be a general trend that the 
radio decays are a bit shallower than the optical / X-ray, so an increasing magnetic energy 
improves the fit. There are other possible causes; any effect that increases the peak flux, 
such as an increasing energy or density, will flatten the radio decay. 

For a sufficiently steep e# oc 7^; x > the resulting model behaviour can no longer fit 

1 11 

the datasets despite any compensating changes in other parameters. As the peak flux oc e B , 
one of the chief spectral behaviour changes is that the peak flux drops, producing steeper 
decays (see Table 2). This leads to a poor radio fit by x ~ +3, often with the radio peak 
flux too low in order to give an appropriate earlier peak flux relative to the optical data. 
The effect is stronger post-jet, and as shown in Table 3 the data for 000926 (with a clear 
jet break in the optical « 2 days) does not even produce a good fit assuming x = +1. The 
extra steepening in the radio requires a later jet break than is compatible with the optical; 
attempts were made to compensate in fits with tj et ^ 3 days, but the radio could not be 
forced to fit. There can also be some fit difficulties at higher frequencies; these may be due 
to a drop in IC flux components for higher early e#, or a change in the spectral slope from 
its link to the decay rates. These effects also combine; a decreased IC flux in the X-ray 
requires an increased synchrotron flux. In some models it is produced by overproducing 
the model's optical flux and suppressing it with host extinction. The combination does not 
produce appropriate spectral indices both in the optical and from the optical to the X-ray. 

While the x — — 1 assumption produced improved fits, for a negative enough 1, oc 7 X 
can no longer match the data, at x ~ —3. For x C 0, will become unphysically high, 
capping at 100%, even with a reasonable early value. This occurs at x = — 4 for 980329, 
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where no good fit can be found with e# < 100% in the data range; early radio flux would 
decline rather than rising, in contradiction to the data (see Table 2). However, the others 
cannot fit before reaching that extreme. By x = —3, 980703's radio peak cannot be matched. 
The peak frequency drops so slowly that it does not pass until quite late; a decline is produced 
post-jet (so too early) by pushing e e — > 1 to maximize the energy losses. This unphysical 
result is also required to counteract the shallow decay at optical and X-ray frequencies. 
970508 fits poorly in a similar manner, but 000926's poor fit is again due to the definite 
optical break. With IC cooling dominant, the pre-jet flux is too shallow, so the jet break is 
pushed back earlier than the data, with a break due to the steepening associated with the 
change to synchrotron cooling dominance which does not match the data well. 

In summary, the data does not constrain the fireball model to have a constant magnetic 
energy fraction e#. Good fits are possible with both increasing and decreasing e#, generally 
as strongly as e B oc 7 +1 through r )~ 2 . With 7 changing by a factor ~ 10 over the data range, 
this allows an extra change in magnetic field strength by a factor of ~3-10. Moreover the 
increasing magnetic energy with e# oc 7 -1 tends to improve the model fit to the radio data. 
As B oc 7e]/ 2 , a more constant level B oc 7 1 / 2 fits slightly better than the basic model. 

5.2. An Investigation of Density Profiles n oc r s 

The CBM density distribution is an important clue to the nature of the progenitor, with 
an r~ 2 profile expected as the signature of a massive star (e.g. Dai & Lu 1998; Chevalier & 
Li 2000). A number of profiles have already been considered. This includes "naked GRBs" 
where a constant density drops to zero after a certain radius (Kumar & Panaitescu 2000) 
and Ramirez- Ruiz et al. (2001 )'s evolution of the afterglow through the ejecta left by a 
Wolf-Rayet star showing that r~ 2 is a very crude approximation to the environment about a 
massive star. Nevertheless, most model fits do as well or better with the ISM approximation 
to the density profile than the Wind one (e.g., Panaitescu & Kumar 2002), although there 
is one case where a Wind CBM fits and an ISM cannot (Price et al. 2002). 

We investigate a wide range of density profiles, parametrized as n oc r s . We recalculated 
the equations for the model presented in §2.1 allowing for a general powerlaw density profile 
index S. This includes the changes to the radiative loss estimates, as in Cohen, Piran & 
Sari (1998). There are no general calculated adjustments in the equation normalizations 
accounting for the postshock electron distribution with a generic density profile, so we used 
the order of magnitude estimates (including for 7, which gives tj et ), and the estimate that 
from the jet expansion tNti/tjet ~ Q 2 ■ 
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We expect cases of S < to help constrain how sensitive the model is to the details 
of a mass-loss wind from the supposed progenitor star. Cases of S < — 2 approach that of 
an evacuated cavity. The equations for the fireball energy from Blandford & McKee (1976) 
break down at S = -4 (see Best & Sari 2000, for S < -4). 

Conversely, S > models a fireball plowing into a medium that gradually increases 
in density. This could mimic a denser region surrounding the burst (though not a sharply 
bounded overdensity) . For S 1, it offers insight into the behaviour when hitting the edge 
of a very dense, but not sharp, shell surrounding the burst. It does not directly model the 
extinction column from the material, so the density would have to cut off before the material 
ahead would absorb all the light emitted from the shock. 

Following is the the expected behaviour of the spectral breaks with general S. v a is for 
the order v a < v m < v c . When v a < v c < v m , it is multiplied by {y m l v c) x ^ 2 \ other 
orderings have other factors. 
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An increasing density (S > 0) does not change the behaviour as greatly as a rarefaction 
(S < 0). With the converse, rapid changes are expected for decreasing S < 0, not increasing 
S > 0. The reasons are illustrated by the range in densities probed and the rate at which 
the shock 7 slows depend upon a density profile. 

Taking the observer time t affected by relativistic beaming, and the relativistic energy, 
E = M'-f 2 (c=l), with M from the density profile n oc r s the results are 

E oc r 5+ V, t oc r^ 2 

Solving for r gives 

n oc r s oc 

As evident from the above equation, for S > a few the rate at which density n increases 
with observer time is weakly dependent upon S, despite its strong dependence on radius. As 
well, in the adiabatic approximation, with E constant 



£-(S+3)/(2S+8) 
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For S > 0, the rate of slowing only varies from from 7 oc t 3//§ to 7 oc t 1 ^ 2 ; it depends quite 
weakly on the density profile for S larger than a few. 

Moreover, in our model, the post-jet shock evolution is not sensitive to density (since 
it is stopped and expands laterally) until the nonrelativistic transition. Changes in the 
parameters to maintain the same jet break for different S may make the nonrelativistic 
transition come at a somewhat different time, but the effect is not dramatic. The question 
becomes whether the model with a given S can reproduce the behaviour pre-jet. To examine 
this, we choose the two datasets with sufficient data for well constrained jet break times. 
GRB 970508, which was well fit with near-isotropy, and GRB 000926, with a sharp jet break 
seen in the optical, are most useful to constrain possible model fits for various S. The others 
had more limited optical data, where jet breaks are generally most obvious, as it was either 
sparsely sampled or partly masked by a bright host. (In §4 we note that the GRB 980703 
dataset could be fit equally well by ISM and Wind-like density profiles; clearly at least some 
very good datasets do not constrain n(r).) 

We have seen cases where n oc r -2 densities can fit the data, but it is difficult to even 
marginally fit models with as steep a gradient as r -2-5 (see the best fits in Table 5). For 
000926, we could not find a satisfactory fit with r -2-5 , as no good fit was possible with a jet 
break. The peak flux declines rapidly prejet; if it scales to fit the early optical it is low by 
the jet break, then too low in the radio. The radio model's decay begins too early due to z/ a 's 
rapid decline as the shock moves into less dense material. GRB 970508 is also not well fit by 
n oc r~ 2 ' 5 . Its r -2 " 5 model with the best x 2 i n Table 5 has an unphysical e# = 100%, giving 
high initial spectral breaks v a and v m as these drop quickly. The peak is high early (as it 
drops quickly); the spectrum from the peak to the optical frequencies is steeper than in the 
ISM case. That steepness out to the X-ray frequencies causes that model to underpredict 
the X-ray flux. The best 970508 model for r~ 2 - 5 under the constraint that e B < 100% does 
not fit the radio data; its predicted flux rises and falls too sharply as u a and v m fall rapidly. 

As mentioned above, for an increasing density gradient the shock slows rapidly compared 
to the constant density case, and the resulting changes in n seen by the shock are gradual. As 
a result, 970508 and 000926 can both be fit with S « 10 (Table 5). For 000926, the increasing 
density leads to a more rapid nonrelativistic transition, improving the fit somewhat relative 
to the constant density case. Further changes will go too far; the fit becomes marginal around 
S — 12 where the best fit models overestimate the early radio flux. 970508 has no visible 
jet break, and is not sensitive to S even for values exceeding S <; 10. Pushing to extreme S 
to discover how far the assumptions can fit the 970508 data (eventually the gradient would 
make the nonrelativistic transition come too early) is of limited value; GRB 970508 's data 
is not very sensitive to an increasing CBM density gradient. 
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In summary, while the data may sometimes accomodate an r~ 2 CBM, it does not fit an 
extreme blown-out density r -2 5 . But a shock plowing into a denser region cannot be easily 
excluded by the data. This is not the same as a sudden jump in density, but a gradual, 
continuous increase n oc r s , for S 1, which is not very realistic but may roughly mimic 
a dense but not sharp shell of material (perhaps ejecta from the progenitor). We conclude 
that the fireball model data fits are not very sensitive to increasing density gradients. 

6. Improving Observational Constraints 

In previous sections (§5.1, §5.2), we found good fits to the data under different assump- 
tions (eg oc 7" 1 , 7 4 " 1 ; n oc r~ 10 , r" 1 ), as reasonable as the good fits derived from the basic 
model of §2.1. The datasets will tolerate widely differing assumptions since the fits need only 
match up with data over a limited range of frequencies and times, and there are significant 
degeneracies between parameters and model assumptions. For example, the decay rate de- 
pends upon the spectral index as well as the assumed n(r) or 65(7). We cross-compared these 
good fits with those of the basic model, checking spectra from 0.01 to 300 days to see the 
best ways of distinguishing them. The fitted models diverge in spectral and temporal regions 
far from the data, as shown in Figures 7 and 8, which highlight the most promising areas 
for improved constraints: longer and more sensitive X-ray and submillimeter observations. 

X-ray lightcurves of some of the acceptable fits are extended to early times for com- 
parison in Figure 7. Sets of acceptable fits for two events are shown to demonstrate that 
the fits in all events diverge by factors of up to three at early times; fluxes of equally ac- 
ceptable fits may be 3 /xJy or 10 /iJy at 0.01-0.03 days. The sensitivity of INTEGRAL'S 
instruments would not be able to distinguish them, but more sensitive X-ray instruments on 
later 7-ray missions such as Swift may. Figure 7 includes a case with an upscattered IC flux 
component, whose peak passage timing gives different curvatures to the light curves. Dense, 
and preferably multifrequency (the IC peak can be seen in the spectrum), X-ray lightcurves 
would break the degeneracy between synchrotron or IC flux as the X-ray source. 

Observations of the broadband peak beyond the radio may be most promising. Figure 
8 shows examples of the submm lightcurves resulting from various acceptable fits found in 
this work; there are a variety of peak levels that subsequently match the radio peak level. 
The peak is in the mid-IR or submm for most of the observable afterglow, from about a 
day to a month. It can rise or fall for a variety of reasons (energy losses, jet break, n(r), 
€5(7)), with details dependent upon the assumptions. These lead to models diverging in 
submm peak height by up to 10 mJy at fractions of a day (~ mJy in the mid-IR). Improved 
submm instruments are expected to reach appropriate sensitivities soon. The ALMA array, 
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to be partially on-line by 2006 and completed by 2010, is expected to give fractional mJy 
sensitivity in a few minutes. Its observations could seriously constrain the peak's behaviour. 

More NIR observations will be of use; in §5.1 we noted that some fits break down where 
host extinction and p could no longer produce appropriate spectral indices both in the optical 
and from the optical to the X-ray. With NIR data, host extinction is better constrained, 
spectral requirements constrain p and allowed model assumptions can be better distinguished 
by their temporal behaviours for that p. 

Finally, earlier optical observations are becoming available now (e.g. GRB 021004 Fox 
(2002), GRB 021211 Fox & Price (2002)) and will be of some use. However, at such early 
times the dominant optical emission should not be due to the synchrotron emission from a 
forward shock into the external medium; reverse shocks (as likely seen in the 990123 optical 
flash, Sari & Piran (1999a); Meszaros & Rees (1999)) and internal shocks may produce 
the early optical flux. This will allow further constraints, but not necessarily to the same 
parameters. We may be seeing in GRB 021004 that the forward shock dominates only after 
ps 0.1 days, and the rise of its peak may be masked by the reverse shock (Kobayashi & Zhang 
2003; Uemura et al. 2003). 

Thus, the behaviours exhibited by good fits (under various model assumptions) to the 
datasets to date may be distinguishable in the near future, especially by densely sampled 
X-ray lightcurves, and observations of the peak at frequencies above the radio. 

7. Conclusions 

We fit four well-studied bursts with extensive radio through X-ray afterglow datasets to 
a fireball model with simple assumptions concerning the microphysics and environment. We 
find a range of reasonable environmental and geometrical parameters. We find all four fit best 
with a constant density medium, one with a value similar to the Milky Way's ISM density, 
n ps 0.2cm" 3 , the other three typical of diffuse clouds n ps 20cm~ 3 . Their kinetic energies 
are comparable to the total GRB 7-ray energy. The collimation varies from near-isotropy to 
a half- angle of 0.04 radians. 

We also find a striking diversity in the fitted microphysical parameter values, far beyond 
the statistical uncertainties. The electron energy distribution index varies from p = 2.1 — 
2.9 and the magnetic energy fraction varies from 0.2% to 25%. As shock physics should 
depend merely on shock strength, we investigated whether the spread could be due to model 
uncertainty, but did not find a set of assumptions which fit the data via universal microphysics 
parameters. 
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We allowed for changes to be made to the model assumptions: eg oc 7 X and indepen- 
dently n oc r s . We find considerable flexibility in the values of x and S that can still produce 
reasonable fits: e# oc 7^; — 2 < x < +1 and n oc r s , with 5 > — 2 through S 1 > 1. More- 
over, some parameter values change by up to an order of magnitude when the assumptions 
underlying the model are altered. Clearly, even the results of very good fits are not unique 
and the parameters cannot be taken at face value. The model assumptions are not strongly 
constrained by the datasets available to date. With this model uncertainty, the evidence 
for massive stellar progenitors from other sources (positions within hosts, possible SN as- 
sociations, see Meszaros 2002) is not hard to reconcile with the lack of clear n oc r~ 2 wind 
signatures in the best fits. Massive stars may not produce a true r~ 2 profile, or its effect 
upon the spectrum could be masked by an inaccuracy in other model assumptions. 

Finally, we compared the spectral evolution of the range of acceptable fits with differing 
assumptions to identify observational strategies that would produce better constraints. Two 
areas are most promising. First, as for now a good fit need only line up with a small 
time range of X-ray observations, the Swift satellite's expected early, well-sampled X-ray 
lightcurves will better constrain the spectral evolution (as well as the IC upscatters of photons 
to the X-ray band and their consistency with the synchrotron model). As well, the peak has 
only been definitively observed at radio frequencies, passing through the mid-IR and submm 
during most of the afterglow. New submm instruments such as ALMA should increase the 
reach of direct peak detections. This will constrain the peak flux evolution, which is sensitive 
to the model assumptions. 

In the future it would be useful to investigate further constraints upon the assumptions. 
These might include variable energy E(t) (under investigation in the context of events such as 
021004, e.g., Heyl & Perna 2003), and the possibility that the electron energy (parameterized 
by e e ), or their acceleration (parameterized by the powerlaw index p), could vary with shock 
strength. 
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Table 1. Fit parameters of the best models with simple assumptions 
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a Time when fast cooling ends at v c — v m 

b Isotropic equivalent blastwave energy (not corrected for collimation), at the time when v c = v m . All tabled energies are in units of 
10 52 ergs, and isotropic-equivalent 

c no lower constraint on this extinction value; 68.3% confidence interval is < 0.037 

Note. — Statistical uncertainties are given for the primary (employed in the fit) parameters; the other columns are derived from the 
fitted values. The quoted uncertainties are produced via the Monte Carlo bootstrap method with 1000 trials to generate the parameter 
distribution. The values bracket the resulting 68.3% confidence interval. These error bars do not include uncertainties in the model 
itself. The model uncertainties are larger than the statistical uncertainties, as demonstrated from the range of parameters that produce 
reasonable fits under various assumptions. 
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Table 2. Model Flux Dependences With es ~ J 



Spectral Region Parameters t, t < tj et a t, t > tj et & 
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a p=2.4 is assumed where necessary 

b for synchrotron cooling dominating; v c behaviour changes for 
IC-dominant cooling, see section 2.2 for details 



Table 3. Fit parameters for the best models with e# = Ky 
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Table 3 — Continued 
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a Time when fast cooling ends at v c = v rn 

isotropic equivalent blastwave energy (not corrected for collimation) , at the time when v c = v m . All 
tabled energies in units of 10 52 ergs, and isotropic-equivalcnt 

c Jet half-opening angle, in radians; > 1 radian treated as isotropic 
d x = quenching rate of energy losses 
c No E(t) quenching, see §5.1 for details 

Note. — Certain items reach limits indicated in the table. Collimations reaching the isotropic limit of 
9 w 1 arc noted as "ISO" . Some magnetic energy fractions reach the physical limit of 100% of the shock 
energy, noted as "MAX"; for such models where the magnetic energy fraction drops with the shock j (x > 0), 
€b is in the physical range over much of the data's time range, the other such models are pinned at this limit 
over the entire dataset. Blank times t cm for the end of fast cooling indicate that the afterglow model formally 
has the transition much earlier than any afterglow data, in the range of prompt emission. 



Table 4. Best fits assuming e B oc 7 
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a Timc when fast cooling ends at v c = u m 

b Isotropic equivalent blastwave energy (not corrected for collimation) , at the time when v c = v m . All tabled energies arc in units of 10 52 ergs, and isotropic- 
cquivalent 

c No results with jet half-opening angles < 1 radian; all treated as isotropic 

d no lower constraint on this extinction value; 68.3% confidence interval is < 0.063 

Note. — Statistical uncertainties are given for the primary (employed in the fit) parameters; the other columns are derived from the fitted values. The quoted 
uncertainties are produced via the Monte Carlo bootstrap method with 1000 trials to generate the parameter distribution. The values bracket the resulting 68.3% 
confidence interval. These error bars do not include uncertainties in the model itself. The model uncertainties are larger than the statistical uncertainties, as 
demonstrated from the range of parameters that produce reasonable fits under various assumptions. 



Table 5. Fit parameters for the best models with n = n^rjr^) 
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a Time when fast cooling ends at v c = v m 

b Isotropic equivalent blastwave energy (not corrected for collimation) , at the time when v c ~v m . All tabled energies are in units of 10 52 ergs, 
and isotropic-equivalcnt 

c Density at R = 10 18 cm, in units of cm~ 3 , the fit parameter 
d in units of cm~ 3 

°Radius in units of 10 18 cm, as calculated, not a fit parameter 
f Jet half-opening angle, in radians 

Note. — Certain items reach limits indicated in the table. Collimations reaching the isotropic limit of 9 w 1 arc noted as "ISO" . In one model 
the magnetic energy fraction reaches the physical limit of 100% of the shock energy, noted as "MAX". The densities given at 1 and 100 days give 
an idea of the range of density probed in the model. These are calculated post-fit; the fit uses scalings based upon the density powerlaw, not a 
density calculation at each time. 



Table 6. Model Flux Dependences With n = n^r/rj) 5 



Spectral Region Parameters t, t < tj et a t, t > t^R 
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a to the model approximation, r is constant during the jet spreading 
phase, and so the behaviour at tj et < t < t NR is the same for any 
density profile (see table 2, for x = 0) 

b for synchrotron cooling dominating; v c behaviour changes for IC- 
dominant cooling, see section 2.2 for details 
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Fig. 1. — Subsets of the optical data from the best fits - the R band light curves for the 
four events. For data selection, see §3. The fits shown are to the full broadband datasets, 
with typically 100+ DOF, and are detailed in Table 1. The GRB 980329 optical data is not 
very constraining for the decay rate - a late optical detection (the early points were found on 
reanalysis) did not allow deep followup over the first week. The scatter in the GRB 970508 
data cannot be explained in any simple model, and otherwise the fits are quite good. 
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Fig. 2. — 8.46 GHz light curves of the 4 events, with the best fits (see Table 1). The light 
grey envelopes show the estimate of the model uncertainty due to interstellar scintillation, 
and data that is not statistically significant at the 2-er are presented as 2-er upper limits 
(downward triangles). As explained in §3 the data for GRB 980329 prior to day four likely 
contains an excess contribution from the reverse shock and was therefore not included in the 
fit. The 980329 radio host component only marginally improves the fit, due to a 1.43 GHz 
average excess, which may indicate a weak radio host flux. The fits are overall quite good, 
although ISS cannot fully account for the scatter in the 970508 dataset (there is, moreover, 
inexplicable scatter at other frequencies). 
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Time Since burst (Days) Time Since burst (Days) 

Fig. 3. — The X-ray data of the 4 events, with the best fits (see Table 1). To parse out 
which model goes with each dataset, the line style is indicated above the symbol labels, and 
linearly-scaled residuals are on the right. The GRB 000926 data is divided into a soft (solids) 
and a hard (open circles with dotted line) band. The broadband fit (with relative flux levels 
and decay rates) indicates an extra flux component in the X-ray, possibly Inverse Compton 
upscatters, but our estimate of this component does not completely fit the data. Better 
sampled X-ray lightcurves, such as those expected from Swift, may clarify deficiencies in 
the X-ray flux model. There is additionally a minor IC flux component in the GRB 980329 
model, and IC dominates the early GRB 980703 X-ray model, providing the flat initial flux 
and slow decay to match the data. 
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Fig. 4. — Parameters from the best fits with the simple model assumptions (see Table 1). 
They are divided into 3 categories: energy & geometry (kinetic energy, collimation half- 
opening angle 9), environment (density) and microphysics (energy partitions: e e for elec- 
trons, €b for magnetic fields, and electron energy distribution index p). See §2 for further 
details concerning the fireball model's parameters. They are presented relative to a nominal 
value as indicated; the error bars are statistical only, 68.3% intervals calculated via Monte 
Carlo bootstraps. The model uncertainties are larger, as seen in the Tables' range of pa- 
rameters that produce reasonable fits under various assumptions. The diversity in energy, 
geometry and environment is not unexpected for some variation in progenitor properties. 
Shock physics, however, is expected to depend only on shock strength. The variation of 
these parameters by orders of magnitude suggests that some effect is unaccounted for in the 
model. 
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Fig. 5. — 8.46 GHz light curves of the four events, with the best fits (solid lines) for an 
assumed magnetic energy relation es oc j x (§5.1). The light grey envelopes are the estimated 
scintillation uncertainties (the 970508 scatter is not fully accounted for, but it has scatter 
excess at all frequencies); data that isn't 3a significant is shown as 2a upper limits (downward 
triangles). 980329 has a radio host component which marginally improves the fit due to a 
1.43 GHz average flux excess. In one case (970508) a constant produces the best fit 
(but only by 1% in the total broadband x 2 ); in the others e# oc 7 _1 gives the best fit. For 
these three, the best constant-e^ fit is shown as a grey dashed line for comparison. The 
model decay for constant-e^ is generally slightly steeper than the data in the late radio; this 
is especially obvious in the last few points of the 000926 light curve. A magnetic energy 
increase at late times flattens this decay, improving the fit. 
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Fig. 6. — Parameters from the best fits assuming e B oc 7 _1 (see §5.1, Table 4). They 
are divided into 3 categories: energy & geometry (kinetic energy, collimation half-opening 
angle 0), environment (density) and microphysics (energy partitions: e e for electrons, e B for 
magnetic fields, and electron energy distribution index p). As it varies, e B is presented both 
with its early (at fast- to slow-cooling transition) and final values. See §2 for further details 
concerning the fireball model's parameters. They are presented relative to a nominal value 
as indicated; the error bars are statistical only, 68.3% intervals calculated via Monte Carlo 
bootstraps. The model uncertainties are larger, as seen in the Tables' range of parameters 
that produce reasonable fits under various assumptions. This model assumption produces fits 
as good as (970508) or better than (980329, 980703, 000926) those assuming e B is constant. 
However, it does not fit with a universality in the microphysics. There is clearly great 
flexibility in the model assumptions allowed by the data, and considerable model uncertainty 
in the derived parameters. 
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Fig. 7. — Comparison of model X-ray lightcurves for several equally acceptable model fits. 
The lightcurves are for a frequency of 6 x 10 17 Hz (nominal for Swift), and show the good fits 
for GRBs 970508 and 000926 with the basic model as well as with the assumption that the 
magnetic energy fraction oc 7 ±1 . While the models are close around the times of the X-ray 
observations (~ days), they diverge at early times, with a spread of several /zJy at 0.01-0.03 
days. Moreover, in the case of 000926, there is significant IC upscattered flux, which gives 
different peak passage times under the differing model assumptions. Early, sensitive, densely 
sampled lightcurves such as those expected from Swift would determine if we are modelling 
the X-ray flux, including the upscattered IC photons, correctly. 
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Fig. 8. — Comparison of model sub millimeter lightcurves for several equally acceptable 
model fits. The lightcurves are for a frequency of 320 GHz (nominal center of an ALMA 
atmospheric window), and show the good fits for GRBs 970508 and 000926 with the basic 
model as well as with the assumption that the density oc r _1 ,r~ 10 . Present sensitivities of ~1 
mJy, attainable only on timescales ~ day, are insufficient to distinguish between the variety 
of peak levels that subsequently match the radio peak. This spread in peak levels is due 
to differing peak behaviours (rising or falling) whose details depend upon factors such as 
energy losses, jet break, n(r), or £5(7). The early model divergences due to density profile 
are of ~ 3 mJy; with models with differing magnetic energy fraction (65(7)) it is up to 
10 mJy. These early differences could be resolved with improved submm instruments soon. 
The ALMA array, to be partially on-line by 2006 and completed by 2010, is expected to give 
fractional mJy sensitivity in a few minutes, which could distinguish amongst these. 



